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We present analytical results for the time-dependent information entropy in exactly solvable two-state (qubit) 
models. The first model describes dephasing (decoherence) in a qubit coupled to a bath of harmonic oscillators. 
The entropy production for this model in the regimes of "complete" and "incomplete" decoherence is discussed. 
As another example, we consider the damped Jaynes-Cummings model describing a spontaneous decay of a 
two-level system into the field vacuum. It is shown that, for all strengths of coupling, the open system passes 
through the mixed state with the maximum information entropy. 
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1. Introduction 

The notion of information entropy has attracted a renewed interest over the last few decades in con- 
nection with fundamental problems in the theory of open quantum systems and quantum entanglement 
(see, e.g., \J,-;^ and references therein). Although the literature concerning diverse aspects of this topic 
is now quite voluminous, not much is known about the entropy behavior in concrete open systems ex- 
hibiting especially intriguing features of quantum dynamics: memory effects, dephasing (decoherence), 
entanglement, etc. 

The simplest models describing many fundamental dynamic properties of open quantum systems 
are two-state systems. Such systems themselves deserve thorough studies as the elementary carriers of 
quantum information (qubits) [5^^. It is also important to note that some of two-state models admit exact 
solutions. The latter fact allows one to gain a valuable insight into general properties of open quantum 
systems. From this point of view, it is of interest to analyze the time behavior of information entropy in 
exactly solvable models. 

In this paper we present exact analytic results for the time-dependent information entropy in two 
physically reasonable qubit models which are frequently used in discussing different problems in the 
theory of open quantum systems. 



2. Entropy of a qubit 

We consider a two-state quantum system (qubit) coupled to a reservoir. In what follows, the qubit 
and the reservoir will be referred to as subsystems A and B, respectively. 

Suppose that at time t the state of the combined system (qubit plus reservoir) is described by some 
density matrix g^^it). Then, the reduced density matrix of the qubit is defined as 

£)^(f) = Tr5£)^B(f), (1) 
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where Tr^ denotes the trace over the reservoir degrees of freedom. The von Neumann (information) 
entropy of the qubit is given by 

S^{t)^-TT^[Q^{t)lng^{t)]. (2) 

It is convenient to use the "spin" representation for a qubit by writing its orthonormal basis |0> and 

|1> as 

I0) = f?l, ll> = W. (3) 



Then, all operators referring to a qubit can be expressed in terms of the Pauli matrices a - {a^,a2,(J^}. 
In particular, the density matrix Q) can be written in the form Rlsl] 

g^{t)^^[l + d-v[t)], (4) 

where 

v{t)^'TT^[ag^{t)] (5) 

is the so-called Bloch vector. In 1 9] we have shown that there exists another representation for the qubit 
density matrix, which is better suited to calculate lng^{t) in (2): 



1 ; ^ — , ^ 1 ll+v 

fAtf) = 2 V 1- (^) exp[a- M(r)], M= - In Y^T^ 



(6) 



Here, v{t) is the modulus of the Bloch vector and uXt v. Strictly speaking, the representation (6) is valid 
only for a mixed state of a qubit with v < 1. Note, however, that the limit v ^ 1 can be taken directly in 
the entropy l|2) after calculating the trace. Using expressions (6), one easily derives from l|2) 

S^(r) = ln2-i(l+y)ln(l+i^)-i(l-i;)ln(l-y). (7) 

For a pure state (y 1), we have - 0, as it should be. The entropy has its maximum - ln2 in the 
mixed state with v-0, when the density matrix is diagonal and (0|p^|0> - (l|p^|l> - 1/2. 
The square modulus of the Bloch vector ([5) can in general be written as 

vHt)^4v^[t)v_[t)+vl[t), (8) 

where 

j^+(f) = Tr^[cr+p^(f)], (9) 

and a+ = [a^± ia^) 12. The two terms in (8) have different physical interpretations. The quantities Vj^{t)- 
<0|p^(f)|l) and v_{t) - <l|p^(f)|0) are often referred to as the coherences. They describe the environ- 
mentally induced dephasing |1]. On the other hand, the time behavior of the component of the Bloch 
vector - (l|p^(f)|l> - (0|p^(f)|0> is determined by energy exchange between an open system and its 
environment, which is responsible for complete statistical equilibrium in the combined system. Thus, for- 
mulas and (8) are convenient for studying the role of different relaxation mechanisms in the entropy 
production. 



3. Entropy production in a dephasing model 

We start with a simple spin-boson model describing a qubit coupled to a reservoir of harmonic oscil- 
lators fll Eoillill . The Hamiltonian of the model is (in our units h-l) 

H^H^ + H^ + Hj^'^a. + Z ^tblbt + '^s E {gA + g*k h] - dO) 

where cJq is the energy difference between the excited state |1> and the ground state |0> of the qubit. 
Bosonic operators fo^ and b^. correspond to the fcth reservoir mode with frequency Wj.. 
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Note that commutes with the Hamiltonian (To). As a consequence, the populations <0|p^(f)|0> 
and (l|p^(f)|l> do not depend on time. In other words, there is no relaxation to a complete equilibrium 
between the qubit and the environment; that is, the model is nonergodic. However, we shall see below that 
this model exhibits a dephasing relaxation and entropy production without energy exchange between the 
qubit and the environment. 

Let us first specify the initial density matrix g^g (0) of the combined system. Usually (see, e.g., (ll 
lioillill ) it is assumed that the subsystems A and B are uncorrelated, and the reservoir is in thermal 
equilibrium at some temperature T. In this paper we will concentrate on the entropy production in time- 
dependent entangled quantum states of the combined system. We assume that at time t-0, the combined 
system is prepared in a pure quantum state which is a direct product 

W^B^O)) = (flfllO) + flilD) «> |0b>, (11) 

where + laj^ = 1, and |0g> denotes the ground state of the reservoir. The initial state 10^) is chosen 
only for simplicity's sake. The subsequent discussion may easily be extended to the case of an arbitrary 
initial state |i//^(0)> of the reservoir. The density matrbc corresponding to (TT) is 

£'ab(0) = IVab(0)><Vab(0)I- (12) 

Since the evolution of the combined system is unitary, the initial state (TT\ evolves after time t into the 
pure state 

W^sit)) = exp{-iHt) ly/^gm), (13) 
so that the density matrix of the combined system is given by 

In principle, the qubit density matrix p^(f), and then the entropy S^(f) can be calculated by using (14). 
It is more convenient, however, to calculate the modulus of the Bloch vector, v(,t), and then apply for- 
mula 0- Since is constant and is determined by the amplitudes and in (11), we need only to 
consider the coherences v^{t). Note that expression (9) can be rewritten as 

v+it]^TT^jj[cF^(t)g^gm], (IS) 

where (T+(f) are the Heisenberg picture operators and the trace is taken over aU degrees of freedom of 
the combined system. In the model l ITO) . equations of motion for a_j_{t) can be solved exactly. The result 
reads \^ 

a+{t)^exp[+icoQt+R[t)](T^, (16) 



where the operator R{t) acts only on the reservoir states and is given by 



r 1 1-e t 
Rit]^j:[^,ml-al{t]b,\, a,{t]^2g,— . (17) 

Substituting the expression (16) into (15) and recalling formula (12) . we find 

y+(f) = v+(0)exp[+iwof-7^a,(f)] (18) 
with the so-called vacuum decoherence function ^ 

r^3,(f) = -ln<05|exp[i?(fl]|05> = -Y,^n{Os\exp[a^{t]bl-al{t]b,^]\0,^). (19) 

k 

After simple algebra which we omit, we obtain 

oo 

/I - COS cot 
dw/(w) , (20) 
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where the continuum limit of the reservoir modes is performed, and the spectral density /(w) is intro- 
duced by the rule 

oo 

E4lgfcl^/(<^fc)= [ do) J {(»)/{(»). (21) 
" 

Now using the solution l ITS) and taking into account that y(0) - 1, we find from equation (8) 

v^it)^ v^ + {l-vl)exp[-2r^^,{t)]. (22) 

Formulas l[7) and (22) determine the time evolution of the qubit entropy. To go beyond these formal 
relations, one needs some information on the spectral density /(w). In many cases of physical interest 
(see, e.g., /(w) may be considered to be a reasonably smooth function which has a power-law 

behavior J{w) oc co^ (s > 0) at frequencies much less than some "cutoff" frequency Q, characteristic of 
the reservoir modes. In the limit co ^ oo, /(&)) is assumed to fall off at least as some negative power of 
0). For the spectral density, we shall take the expression which is most commonly used in the theory of 
spin-boson systems U(Vil3i1 : 

/M = A.O^-'w'e"'"'", (23) 

where Aj is a dimensionless coupling constant. The case 5 = 1 is usually called the "Ohmic" case, the case 
5 > 1 is "super-Ohmic", and the case < 5 < 1 is "sub-Ohmic". 

Substituting the spectral density (23) into l [20) and doing standard integrals, one gets 



, „, ,,L cos[(5-l)arctan(Or)] 1 
A 

rvac(f) = y In(l + a2f2), (5=1), 

where T(s) is the Euler gamma function. Some important properties of 7vac(^) can easily be seen directly 
from the above expressions. First, Jvac ( f ) is a monotonously increasing function of time for 5^1. Second, 
in the super-Ohmic case 7vac(fl has a long-time limit: 

7^3e(oo)= limrvacU) = A,r(s-l), (s>l). (25) 

oo 

Finally, the 7vac(^) monotonously saturates to Jvac^oo) for 1 < s ^ 2 and is a nonmonotonous function of 
time for s > 2. 

In discussing the properties of the qubit entropy in the model tlO) . it is necessary to distinguish two 
cases: the regime of "complete decoherence", and the regime of "incomplete decoherence". In the former 
case (s ^ 1) we have 7vac(^) ^ oo as f ^ cxd, and hence (f) ^ 0. Then, it follows directly from l[7) and l l22t 
that the limiting value of the qubit entropy is given by 

S^(oo) = ln2-i(l-Hy3l)ln(l-Hk3l)-i(l-|i;3l)ln(l-|i;3l), [s^l). (26) 

The maximum qubit entropy S^^ax^"") = ln2 corresponds to the initial state with equal populations (v^ - 
0). In the case of "incomplete decoherence", we have 

S^(cx,) = ln2-i(l+i;Jln(l+i;J-i(l-yJln(l-yJ, (5>1), (27) 

where 

Voo^{vl + {l- vl) exp [-27^,, {oo)]f\ (28) 

Figure[l]illustrates the time behavior of the qubit entropy for different values of the parameter 5 and for 
different coupling strengths. It should be emphasized at once that, for the case under consideration here, 
the entropy production has no "thermodynamic" meaning. Indeed, at any time f, the combined system is 
in the pure quantum state (13). Note, however, that this state is entangled |1], i.e., it cannot be written as 
a direct product of states of the subsystems. Thus, the information entropy S^(f) may be regarded as a 
measure of entanglement. 
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Figure 1. Time evolution of the qubit entropy in the dephasing model (W\ . In all cases A = A^. The qubit 
is initially in the pure state with equal populations {v^ = 0). 



Within the framework of the model (TO) with the initial condition (TD . the structure of the entan- 
gled state (T3\ can easily be established. Using the fact that the Hamiltonian (lO) contains only the qubit 
operator 03, we obtain 

ly/^eit)) = floe-Hf'2|^H(f)^ ^ + ^^gH*/2|^W(j)) ^ n) (29) 

with 

\yf^^\t))^exp{-iH^^h)\0^), H^B-^^L'^Ah^Llskbl + glbX (30) 

k k 

The reservoir states |i/'^'(f)> are normalized but {y/^j^Ht)\y/^j^\t)) + 0. It is possible, however, to write the 
states appearing in l [29) as linear combinations (the argument t is omitted for brevity) 

IVb^> = ;6o"^I*o> + M"'|1'i>' |O> = ao|(/)o> + ail0i>, |1> = a[,|(/)o> + (31) 

where the amplitudes may be determined from the following conditions: a) the pairs of states (|<5o)'l'l'i)) 
and(|(/)Q>,|(/)^>) are orthonormal; b) the terms with \<b^®\(p^ and |<1>j>«i|0q> cancel when expressions 
are inserted into l[29). Then, the state of the combined system takes the form 

|i/'^B(f)> = Ao(f)|<Do(f)> «. |(/)o(f)> + ^i(f)|(l)i(f)> % |(/)i(f)> (32) 

with amplitudes satisfying \AQ{t)\^ ^\A^{t)'^ -\. Explicit expressions for Ag(f) and A^{t) are rather cum- 
bersome and are not given here. 

Formula (32) is an example of the so-called Schmidt decomposition of quantum states of a combined 
system Amplitudes Ag(f) and A^[t) play a role of the corresponding Schmidt coefficients. From (32) it 
follows that the reduced density matrices of the subsystems may be written as 

gA{t)^\Ao\^\(Po)(cPo\ + \A^\^\cp^)(cp^\, gB(t)^\Af\%){%\ + \Af\<^^)(<t>^\. (33) 

Using these formulas, it is easy to show that 

S^it) = Sgit) = -\AQ[t)\Hn\Ao[t)\^-\A^[t)\Hn\A^{t)\\ (34) 

which is the well-known consequence of the Schmidt decomposition theorem for entangled quantum 
states ill]. 



4. Spontaneous decay of a qubit 

Now we consider the entropy production in an exactly solvable model which describes a spontaneous 
decay of a two-level system into a field vacuum (iIQEID- The total Hamiltonian of the model is 

H^H^ + H, + H,^^a, + Y.'^Ah + L [Sk^l^- + Sl\cT,] , (35) 
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where the index k labels the photon modes with frequencies co^. 

The initial state of the combined system is again assumed to be given by JTT) . i.e., i;(0) - 1 and, conse- 
quently, S^(0) - 0. To obtain the entropy of the qubit at times t>Owe need to calculate y(f). To do this, 
we will closely follow the approach taken in the works 

It is convenient to work in the interaction picture with the unperturbed Hamiltonian - Hj^ + H^. 
Introducing the interaction picture state vector |i//^g(r)) - exp(iHQr) |i/'^g(?)>, and applying the standard 
procedure, one obtains 

t 



IVab(^)> = exp+ -i j dt'Hjit') 



where exp+ [■ ■ ■ ] is the chronologically ordered exponent, and 



IV^b(0)>, (36) 



k 

is the interaction picture Hamiltonian. As noted in the work (see also fl'.'li]), there exists a simple 
representation for the state (36), which can be derived using the following properties of the Hamilto- 
nian (37): 

Hjit) |0> » lOjj) = 0, H^(r) |1> » \0s) = Es'fce-'"""""*'* |0> » ^lO^}, 

k (38) 
H,{t] |0> » bl\Os) = g*e't'"»-'-^-'' |1> ^ lO^}. 

Recalling the expression dlT) . one can easily show that at any time f, the state l[36) is a superposition of 

10) «. lO^), |1> «. |0g), and |0> ® fo|.|Ojj): 

IV^^B(f)> = [flfllO) + q (DID] «. \0g) + ^ Cj.(r)|0> <S> fo[jOB>. (39) 

The amplitudes Cj(t) and Cj.(r) satisfy the initial conditions Cj(0) = a^, Cj.(0) = 0, and the normaUzation 
condition 

\af + \c^{t]\^ + Y,\c,,it)\^^l. (40) 

k 

The qubit density matrix in the interaction picture can now be calculated using l [39) : 

( f) = Tr^ { 1 ( f) > ( f) I }= i + i (2 1 Ci ( f) 1 2 - 1) cTg + Ci ( f) cr+ + flo Ci ( f) cr_ , (41) 



where the amplitudes Cj.(f) have been eliminated with the help of 00). Expressions for v_j_{t) and Vjtf) 
follow immediately from l[5) and the relation g^{t) - exp{-iHAt)g ^(t) exp(i/fyif): 

i^+(f) = e''^o'floq*(r), y_(f) = e"'"o*a*q(f), i^glf) = 2|Ci (f)!^ - 1. (42) 

Substituting these expressions into l|8) and using la^f = 1 - |Cj (0)|^, we finally obtain 

i;2(f)= l-4|Ci(f)|2[|ci(0)|2-|q(f)|2]. (43) 

Thus, the modulus of the Bloch vector and, consequently, the qubit entropy S^(f) are completely deter- 
mined by |Cj(f)|^ which is the time-dependent population of the excited state |1). It is important to note 
that the probability amplitude c^[t) obeys a closed equation fllEill 



T 

c^[t)^-j At' f(t-t')c^{t') 



k J 



(44) 



with the kernel 



43004-6 



Entropy production in open quantum systems: exactly solvable qubit models 



where /(w) is the spectral density of the field. In some important special cases of /(w), equation (44) can 
be solved to give exact analytic solutions for c^it). We shall restrict ourselves to the so-called damped 
Jaynes-Cummings model (see, e.g., (llllil]) which describes the resonant coupling of a two-level atom to a 
single cavity mode of the field. In this model, the effective spectral density has the form 



2^ (wq-w) +A2 



(46) 



where A is a spectral width of the coupling, and the parameter 7q defines the characteristic time scale 
= I/Tq on which the state of the qubit changes. For the details of solving the equation j44) we refer 

toiillii.i 



Here, we quote the resulting expression for the population ICj^tf)!^: 



|ci(f)|2 = |ci(0)|2e-^* 



A 1^ 
cosh(Ar/2)-H — sinh(Af/2) , 
A 



(47) 



where A = (A^ - Zy^A)^'^. Substituting the above expression into (43) and introducing the dimensionless 
coupling parameter 



we obtain 



vHt) = l-4|Ci(0)|*e-'^^F^(r) [l-e-^'Fjf(r)] . 



(48) 



(49) 



The form of the function PAt) depends on the value of K: 



FAt) ■■ 



cosh i Vl-KAt/z] + ^ sinh ( Vl-KAtl2\ 
cos [ ^/K-\ltl2)^ + sin [ VK-lAt/2^ 



K< 1, 

K> 1, 
1. 



(SO) 



Formulas and (49) completely determine the qubit entropy in the Jaynes-Cummings model. Figure [2] 




b) 




Figure 2. Time evolution of the qubit entropy in the damped Jaynes-Cummings model: a) Weak and mod- 
erate coupling {K < 1); b) Strong coupling (K > 1). In all cases the qubit is initially in the excited state |1>, 
i.e., |Ci(0)|2 = l. 



illustrates the time behavior of S^(f). We take the situation where the initial pure state |1> of the qubit 
evolves into the final pure state |0). Note that for all strengths of coupling, at some time t^^, the qubit 
passes through the mixed state with the maximum entropy S^{t^) - in 2, i.e., with v{tj^) - 0. It is inter- 
esting that this "maximally entangled" intermediate state of the combined system appears only in the 
case of the maximum initial population of the excited qubit state 1 1> (| Cj (0) 1^ = 1), as may be seen directly 
from (43). The time behavior of the qubit entropy for smaUer values of |Cj(0)|2 is shown in figure [S] 
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Figure 3. Time evolution of the qubit entropy in the damped Jaynes-Cummings model for different initial 
populations of the excited qubit state 1 1> in the cases of weak and strong coupling. 

5. Conclusion 

In this paper we have considered the time evolution of information entropy in two exactly solvable 
models of two-state open quantum systems (qubits). Calculations were based on the simple but quit gen- 
eral representation Q for the qubit entropy. 

Our discussion was restricted to the special case of the qubit and the reservoir being initially uncor- 
related and the reservoir being in the ground state, i.e., at zero temperature. This simplifying assumption 
allows one to study in detail the entropy production in entangled time-dependent quantum states of a 
combined system. It should be noted, however, that in many situations of physical interest, a factorized 
(uncorrelated) initial state at T = cannot always be realized, so that the dynamics of thermal and cor- 
related initial states, including the entropy behavior, is of great significance. For instance, it was shown 
in 1 9] that the dephasing model (TOt admits exact solutions for a large class of physically reasonable cor- 
related initial states at finite temperatures. It was found that, for a sufficiently strong coupling, initial 
qubit-environment correlations have a profound effect on the dephasing process and on the time be- 
havior of entropy. It would be of interest to study the entropy behavior in further examples of exactly 
solvable models of open quantum systems in the presence of initial system-environment correlations. 
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ripoAyKyBaHHfl eHTponii y bIakpi/iti/ix KBaHToei/ix ci/icreiviax: 
TOMHO po3B'5i3Hi Ky6iT MO/\en\ 

B.r. Mopo30#], r. PenKeP 

1 ^ ^ ^ 

MocKOBCbKnn flepxaBHwn TexHHHwn yHiBepcMTer pafliorexHiKu, e/ieKipoHiKu la aBTOMaii/iKu, 
npocneKT BepHaflCbKoro, 78, 119454 MocKBa, Pocia 
^ YHiBepcuTeT m. Poctok, YHiBepcuTeTCbKa n/iou^a, 3, D-18051 m. Poctok, HiweHHUHa 

Mm npeflCTaB/i^ieMO aHa;iiTH4Hi pesy/ibTaTH flyia 4acoB03a;iexHOi iH(()opMai4iMHoT eHiponii' b flBOcraHOBHx to- 
HHO po3B'fl3Hnx {Ky6iT) MOfle/iax. flepuja MOfle/ib onucye flecfiasyBaHH^ (fleKorepeHL4iK)) b Ky6iTi, akhm e sb'a- 
saHni^ 3 pe3epByapoM rapMOHi^Hnx 0CL4M;iflT0piB. 06roBopK)eTbCfl npoflyKyBaHHfl eHiponiT fl/ia i^iei MOfle/ii y 
pexuMax "noBHoT" ra "HenoBHoT" fleKorepeHi^ii. Ihujum npuK/iafl mm po3r;iflflaeMO 3afleMnc|)0BaHy MOfle/ib 
flxeMHca-KaMMiHrca, ^Ka onucye cnoHiaHHe 3aTyxaHH5i flBoplBHeBoT CMCreMU b no/ibOBOMy BaKyyMi. noKa3a- 
HO, w,o p,n!f Bcix cvin 3B'fl3Ky, BlflKpnTa CMcreMa nepexoflHTb 4epe3 3MiujaHMM craH 3 MaKCHMyMOM iH(|3opMaL4lM- 
HOT eHTponii. 

K/iioMOBi c;ioBa: iH<l>opMau,ii/iHa enTponip, BiflKpmi KBaHTOBi cucreMn, Ky6iT MOfl,eni, flBKorepeHLiifi, KsaHTOse 
sanjiyrysaHHfi 
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